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Abstract 

Previous work on t — b — r Yukawa-unified supersymmetry, as expected from SUSY 
GUT theories based on the gauge group 50(10), tended to have exceedingly large 
electroweak fine-tuning (EWFT). Here, we examine supersymmetric models where 
we simultaneously require low EWFT ("natural SUSY") and a high degree of 
Yukawa coupling unification, along with a light Higgs scalar with ~ 125 GeV. 
As Yukawa unification requires large tan/3 ~ 50, while EWFT requires rather light 
third generation squarks and low [i ~ 100 — 250 GeV, 5-physics constraints from 
BR(i? — > X s ^f) and BR(i? s — > /i + /x~) can be severe. We are able to find models 
with EWFT A < 50 - 100 (better than 1-2% EWFT) and with Yukawa unification 
as low as i? yu k ~ 1.3 (30% unification) if .B-physics constraints are imposed. This 
may be improved to -R yu k ~ 1-2 if additional small flavor violating terms conspire 
to improve accord with S-constraints. We present several Yukawa-unified natural 
SUSY (YUNS) benchmark points. LHC searches will be able to access gluinos in the 
lower 1 — 2 TeV portion of their predicted mass range although much of YUNS pa- 
rameter space may lie beyond LHC14 reach. If heavy Higgs bosons can be accessed 
at a high rate, then the rare H, A — >• u + u~ decay might allow a determination of 
tan (3 ~ 50 as predicted by YUNS models. Finally, the predicted light higgsinos 
should be accessible to a linear e + e~ collider with yfs ~ 0.5 TeV. 

Keywords: Supersymmetry Phenomenology, Supersymmetric Standard Model, Large 
Hadron Collider 
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1 Introduction 



A striking feature in nature is that all the fermions of each generation fill out a complete 
16- dimensional spinor multiplet of the gauge group SO (10) pQ. While ordinary grand 
unified theories (GUTs) suffer from the notorious gauge hierarchy problem, supersym- 
metric (SUSY) GUTs not only tame this hierarchy problem [2j, but they also receive 
support from the well-known unification of gauge couplings [3]. In the simplest SO (10) 
SUSY GUT theories, where both MSSM Higgs doublets H u and H d occupy the same 
10-dimensional representation, one also expects unification of third generation Yukawa 
couplings ft, fb and f T at M G ut - 2x 10 16 GeV [U [5]. The t — b — r Yukawa coupling 
unification is highly sensitive to both 2-loop renormalization group running (RGEs) and 
to threshold corrections when transitioning between MSSM and SM effective theories at 
the SUSY particle mass scale. Thus, the entire SUSY mass spectrum enters into a precise 
computation of Yukawa coupling unification. 

Many groups have explored t — b — r Yukawa unification (YU) in SUSY theories [EJCFJIEl 
El [TOl El [121 [TJ]. It has been found that, for /i > 0, YU can occur at the few percent level 
in either the Higgs splitting (HS) model or the DR3 model (D-term splitting, right-hand 
neutrino effects and third generation splitting), provided that the GUT scale soft SUSY 
breaking (SSB) terms are related as 

Al ~ 2mj ~ 4m 2 6 , (1) 

with Aq < 0. Moreover, the GUT scale Higgs splitting m 2 H < m 2 H is needed to allow 
for an appropriate radiative breakdown of electroweak symmetry. In these models, first 
generation squarks and sleptons are required to be in the multi-TeV range [8j [10] while 
third generation sfermions are driven to much lighter TeV-scale masses. A benefit of 
these models is that the light SUSY Higgs boson mass rrih tends naturally to be in the 
125 GeV rangeQas required by the recent LHC discovery jX6|, [T7j. The gaugino masses 
from YU SUSY are expected to be quite light, with typically m g < 500 GeV (now ex- 
cluded by LHC searches for gluino pair production along with cascade decay into states 
containing 6-quarks[IH]), although solutions with heavier gluinos ~ 1 — 2 TeV can also 
be found [TS]. In all these cases, the superpotential /x parameter, which is extracted from 
the electroweak minimization conditions, occupies values in the 1-10 TeV regime, leading 
to severe electroweak fine-tuning (EWFT). 

In this work, we examine to what extent it is possible to reconcile YU with low EWFT. 
Minimization of the SUSY scalar potential allows one to relate the Z mass scale to the 
superpartner mass scale via the well-known relation 

1 M 2 _ ( m H d + Z d ) - (m 2 Hu + S„) tan 2 g 2 

2 Z ~ (tan 2 - 1) ■ {2) 

The radiative corrections H u and are given in the 1-loop approximation of the Higgs 
effective potential by: 

v _ 1 dAV 

v u a on u d 



1 The large negative A$ in Eq. ^ leads to maximal stop mixing, see e.g. OQ2], thus increasing rrih 
to the desired range. 
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where AV is the one-loop correction to the tree-level potential, and the derivative is 
evaluated in the physical vacuum: i.e. the fields are set to their vacuum expectation 
values after evaluating the derivative. At the one-loop level and in the limit of setting 
first/second generation Yukawa couplings to zero, S u contains 18 and contains 19 
separate contributions from various particles/sparticles [20J. We include contributions 
from W ± , Z, t lj2 , &i,2, f i,2) Wi,2, -^1,2,3,4) ti b and r, h, H and H ± . We adopt a scale 
choice Q 2 = m^m^ to minimize the largest of the logarithms. The dominant contribution 
to the terms E M d arise from superpotential Yukawa interactions of third generation squarks 
involving the top quark Yukawa coupling. For instance, the dominant contribution to S„ 
is given by 



3 „, a 



£ u (i 1>2 ) = ^« 



, 2 2 f?A 2 t - 8g 2 z (j - lx w )A t 
It - 9z+ Zl2 



m~ — m~ 



(4) 



where A t = (m~ — m 2 ^)/2 + m| cos2/3(| — \x w ), g\ = {g 2 + g' 2 )/S and x w = sin 2 6 W and 

F{m 2 ) = m 2 (log ^ — 1). This expression thus grows quadratically with the stop mass. 

We adopt the fine-tuning measure from [20], which requires that each of the 40 terms 
on the right-hand-side (RHS) of Eq. pi) should be of order ~ m 2 z /2. Labeling each term as 



d (with i = H d , H u , n, etc.), we may require C max = max|Cj| < A^, 

where A max ~ 100 — 300 GeV, depending on how much EWFT one is willing to tolerate. 
This measure of fine-tuning is similar to (but not exactly the same as) Kitano-Nomura [21] 
but different from Barbieri-Giudice [22] beyond the tree-level. In the following, we will 
use the fine-tuning parameter 

A = C max /(m 2 z /2), (5) 

where lower values of A correspond to less fine-tuning, and e.g. A — 20 would correspond 
to A -1 = 5% fine-tuning. 

Our goal in this paper is to search for parameter choices which 

1. Maximize the degree of Yukawa coupling unification, i.e. minimize 

_ max(/ ; , f b , f T ) 
^ yuk min(/ t , f b , f r ) [b) 

with each Yukawa coupling evaluated at the GUT scale. Thus, a value of i? yu k = 1 
would give perfect Yukawa coupling unification. 

2. Have as low EWFT as possible (in practice, we will require A < 100, or better than 
1% EWFT). 

3. Have rrih ~ 125 GeV in accord with the recent LHC discovery of a Higgs-like reso- 
nance. In practice, we will require 122 GeV < rrih < 128 GeV to allow for a roughly 
2-3 GeV error in the RG-improved one-loop effective potential calculation of the 
Higgs mass rrih. 

We recognize that in addition to EWFT, there also exists a fine-tuning associated with 
generating particular weak scale SUSY spectra from distinct GUT scale parameters [22] , 
(see also [23l El] for related discussions). Here we adopt the less restrictive weak scale 
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fine-tuning condition, which nonetheless turns out to be indeed very restrictive. In this 
vein, we regard particular GUT scale parameters as merely a parametrization of our 
ignorance of the mechanism of SUSY breaking and soft term generation]^] For instance, 
in this paper we will require rather low values of superpotential /x parameter to avoid 
excessive EWFT in Eq. Q. In generic SUSY models, the value of /i is expected to be of 
order Mpi anc k since it is a dimensionful SUSY-preserving parameter. Excessively large \x 
can be avoided ala Giudice-Masiero [26] where the /i superpotential term is forbidden by 
some high scale symmetry, but then is regenerated as a soft SUSY breaking term via a 
Higgs-Higgs coupling to the hidden sector. 

For the remainder of this paper, in Section [2] we present details of our scan over SUSY 
parameter space, and which constraints are invoked in our analysis. In Section [3] we 
present the results of our parameter space scans. We will find that requiring A < 100 and 
rrih ~ 125 GeV only allows for -R yu k as low as ~ 1.2 — 1.3. While this degree of Yukawa 
unification is not optimal, we feel it is still useful in that it might guide model builders 
towards models including additional GUT scale threshold corrections or extra matter or 
above- GUT-scale running which may ameliorate the situation. In Section [5j we discuss 
observable consequences of YUNS for LHC, ILC and dark matter searches. We pay some 
attention to methods which might allow one to distinguish YUNS from generic NS models 
at lower tan /3 values. In Section [6] we present a summary and conclusions. 

2 Parameter space and Yukawa unification 

For our calculations, we adopt the IS A JET 7.83 [27] SUSY spectrum generator ISAS- 
UGRA [2H]- ISASUGRA begins the calculation of the sparticle mass spectrum with input 
DR gauge couplings and /&, f T Yukawa couplings at the scale Q = Mz (ft running be- 
gins at Q = m t ) and evolves the 6 couplings up in energy to scale Q = Mgut (defined 
as the value Q where g\ = #2) using two-loop RGEs. At Q = Mgut, we input the 
soft SUSY breaking parameters as boundary conditions, and evolve the set of 26 cou- 
pled MSSM RGEs [2H] back down in scale to Q = M z . Full two-loop MSSM RGEs are 
used for soft term evolution, while the gauge and Yukawa coupling evolution includes 
threshold effects in the one-loop beta-functions, so the gauge and Yukawa couplings tran- 
sition smoothly from the MSSM to SM effective theories as different mass thresholds are 
passed. In ISASUGRA, the values of SSB terms of sparticles which mix are frozen out at 
the scale Q = Msusy — yJ m t L m t R 5 while non-mixing SSB terms are frozen out at their 
own mass scale [28J. The scalar potential is minimized using the RG-improved one-loop 
MSSM effective potential evaluated at an optimized scale Q = Msusy which accounts for 
leading two-loop effects [30]. Once the tree- level sparticle mass spectrum is computed, 
full one-loop radiative corrections are calculated for all sparticle and Higgs boson masses, 
including complete one-loop weak scale threshold corrections for the top, bottom and tau 
masses at scale Q = M SUSY [21] • Since the GUT scale Yukawa couplings are modified 
by the threshold corrections, the ISAJET RGE solution must be imposed iteratively with 
successive up-down running until a convergent sparticle mass solution is found. Since 
ISASUGRA uses a "tower of effective theories" approach to RG evolution, we expect a 
more accurate evaluation of the sparticle mass spectrum for models with split spectra (this 

2 The relation between GUT scale fine-tuning and the b-r Yukawa coupling ratio was studied in [25] . 
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procedure sums the logarithms of potentially large ratios of sparticle masses) than with 
programs which make an all-at-once transition from the MSSM to SM effective theories. 
The fine-tuning measure A described in Sec. [ljhas been implemented in ISAJET 7.83 [27] . 

In models of "natural SUSY" (NS) |321|2IIE31E11E3ISIEZI|2Q1, the first requirement 
to gain a low EWFT is that the /i parameter be of the order of ~ Mz, while in Yukawa- 
unified SUSY, it is necessary to invoke some manner of Higgs soft term splitting at the 
GUT scale in order to obtain raditive EWSB. For these reasons, the model parameter 
space chosen in this study is the two-parameter non-universal Higgs model (NUHM2) 
where the weak scale values of \i and tua are input in lieu of the GUT scale values of m 2 H 
and m 2 H . In addition, to allow for (sub)TeV-scale third generation masses (as required 
by EWFT from Eq. ^ along with at least a partial decoupling solution to the SUSY 
flavor and CP problems, we allow for split first /second and third generations at the GUT 
scale. Thus, the parameter space we choose is given by 

m 16 (l,2), m 16 (3), m 1/2 , A , tan/3, //, m A . (7) 

Here, m 16 (l, 2) and m 16 (3) are the first/second and third generation sfermion soft masses, 
respectively; = M\ = M 2 = M 3 is the universal gaugino mass parameter; and 

A Q = A t = Af, = A T is the universal trilinear coupling. These parameters are defined at 
Mqut, while tan/3, /x, and m A are defined at the weak scale. The top quark mass is set 
to mt = 173.2 GeV. 

We search for mass spectra with low EWFT A and low -R yu k by performing a vast 
random scan over the following parameter ranges (masses in GeV units): 

5000 GeV < m 16 (l,2) < 20000 GeV, 
< mi 6 (3) < 20000 GeV, 
300 GeV < mi/j < 2000 GeV, 

-2 < A)/m 16 (3) < 1, (8) 
50 < tan/3 < 60, 
100 GeV < /I < 250 GeV, 
500 GeV <m A < 5000 GeV. 

The lower limit on m^j comes from the approximate LHC bound of irtg > 900 GeV (for 
rriq ^> rrig) [SB], while the lower bound on tua comes from LHC searches for A, H — > t + t~ 
which require m A > 500 GeV at tan /3 ~ 50 • We require of our solutions that 

1. electroweak symmetry be radiatively broken (REWSB), 

2. the neutralino Z\ is the lightest MSSM particle, 

3. the light chargino mass obeys the rather model independent LEP2 limit that > 
103.5 GeV HQ], and 

4. the light Higgs mass falls within the window = 122-128 GeV, where we adopt 
±3 GeV as theoretical error on the Higgs mass calculation. 



4 



1.5 
1.45 
1.4 

^ 1.35 
>- 

= 1.3 
1.25 



1.2 



1.15 

100 200 300 400 500 600 700 800 900 1000 

A 

Figure 1: Scatter plot of -R yu k versus A from the scan defined in Eq. (2.2) The dark blue 
triangles violate the i?-physics constraints of Eq. (2.3) The light blue triangles obey the 
BR(B S — > /i + /U~) constraint, but deviate from the measured BR(B — > X s 7) by more 
than 3a. Finally, the pink squares satisfy both the BR(£> S — > //+ and BR(B — > X s ^) 
constraints. Only points with A < 1000 and R yu k < 1.5 are shown. 

Regarding 5-physics constraints, we consider BR(B — > X s ^) = (3.55 ± 0.34) x 
10~ 4 , where experimental and theoretical uncertainties have been added in squares, and 
BR(B S < 4.2 x 10~ 9 at 95% CL. In the following we will impose the BR(B ->• 

X s ^f) constraint only at the 3a level (the reason for this is that calculating BR(B — > X s ~/) 
with both IsaTools and SuperlSO, we observe deviations of the order of 5% in the relevant 
region of parameter space). For the BR(B S — > constraint, we assume a theoretical 

uncertainty of 20%. This leads to the the following limits 

BR(B X s7 ) = [2.53, 4.57] x 10~ 4 , 

BR(B S //+//-) < 5.04 x 10~ 9 , (9) 

which we will use throughout the numerical analysis. 

Regarding neutralino relic density, we remark here that models of natural SUSY con- 
tain a higgsino-like lightest neutralino with thermal abundance of typically Q% p h 2 ~ 
O(10~ 3 — 10~ 2 ). This thermal under-abundance can be regarded as a positive feature 
of NS models in the sense that if one invokes the axion solution to the strong CP prob- 
lem, then one expects mixed axion-higgsino dark matter, where the higgsino portion is 
typically enhanced by thermal axino production and decay to higgsinos in the early uni- 
verse. Thus, a thermal under-abundance leaves room for additional non-thermal higgsino 
production plus an axion component to the dark matter [41]. 

3 Scan results 

As our first result, we show in Fig. [Tjpoints from our parameter space scan in the -R yu k vs. A 
plane. All points have = 122-128 GeV and obey the current LEP and LHC SUSY 
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Figure 2: Dependence of -R yu k on my^, mi 6 (3), A /mi 6 (3) and tan/3, for A < 100. Same 
color code as in Fig. [I] 



mass limits. The dark blue triangles however violate the 5-physics constraints of Eq. (|9]). 
These points may still be valid if additional small flavor violating terms are allowed in 
the theory. The light blue triangles obey the BR(£> S — > constraint, but deviate 

from the measured BK(B — > X s j) by more than 3a. Finally, the pink squares satisfy 
both the BK(B S — > and BR(I? — > X s / y) constraints. This color scheme is used 

throughout the remainder of the paper. We see already from this plot that flavor physics 
constraints significantly affect the parameter space of Yukawa-unified natural SUSY. This 
is to be expected, since naturalness requires lighter third generation squarks while Yukawa 
unification requires tan/3 ~ 50: both these effects bolster SUSY contributions to B- 
physics observables. (Analogous observations were made in [10J H2] in the context of 
generic YU models.) Without flavor constrains, we can obtain i? yu k as low as ~ 1.2. 
The BR(£> S — > /U + /i~) constraint pushes this up to R yu k > 1.27, and the BR(£> — > X s ^) 
constraint to R yu k > 1.3. Aside from the .B-physics constraints, it is intriguing that points 
with lowest A also have lowest values of i? yu k- This is because low A requires light third 
generation squark masses, while at the same time Yukawa unification requires large SUSY 
threshold corrections which also require lighter third generation squarks. 

In Fig. |2j we show the value of -R yu k versus various input parameters for A < 100 (less 
than 1% fine-tuning). Here, the behavior deviates considerably from t — b — r unified 
models with large \x > and a spectrum derived from the radiatively driven inverted 
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scalar mass hierarchy [7J El OH]- From frame a), we see that YUNS models actually 
prefer large mi/ 2 whereas generic Yukawa unified SUSY models (with -R yu k < 1-1 but 
arbitrary EWFT), YUS for short, prefer low m\/2- The preference for large mi/2 helps 
to avoid LHC constraints on the gluino mass. If we impose the I?-physics constraints, 
then the distribution flattens out with some preference for lower mi/2 values. In frame b), 
where -R yu k is plotted vs. m 16 (3), we see that lowest -R yu k values prefer mi 6 (3) ~ 2 TeV, 
which leads to rather light third generation squarks and typically violation of I?-physics 
constraints. If we respect 5-constraints, then larger values of mi6(3) ~ 3 — 7 TeV are 
preferred, at the cost of larger values of -R yu k- In frame c), we plot versus A /mi 6 (3). For 
YUS models, there is a strong preference for Aq ~ — 1m\§ [7J, while for YUNS, the lowest 
-Ryuk values are obtained for smaller \A§\. Imposing S-constraints, the preference moves 
tO A rsj (-2 to 0) x m,xg(3). Frame d) shows i? yu k vs. tan/3. As expected, tan/3 ~ 50 is 
preferred. 

The dependence of A on the input parameters is shown in Fig. [3] for points with 
i? yu k < 1.4. In frame a), we see a mild preference by A for low m^, i.e. the gluino 
mass can't be too heavy, lest it pushes the stop masses too high, leading to large E u . 
From frame b), we see low A prefers mig(3) ~ 2 — 4 TeV. If mig(3) is much higher, 
then third generation squarks are too heavy to give low EWFT, while if mi6(3) is too 
ight, we generate tachyonic spectra: the optimal corresponds to mi6(3) ~ 2 — 4 TeV. For 
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Figure 4: Dependence of -R yu k on rrih, Trig, m~ tx and uia, for A < 100. Same color code as 
in Fig. [Tj 



mi6(3) ~ 1 — 2 TeV, then sub-TeV top squarks are generated leading to violation of B- 
constraints. In frame c), we see that low A allows a wide range of Aq unless S-constraints 
are imposed, in which case A ~ (—2 to 0) x mi 6 (3) is again preferred. The large negative 
Aq values lead to larger m/, values and also can lower the EWFT [23]. In frame d), we 
see that there is only mild preference for tan (3 ~ 48 — 58 values unless S-constraints are 
respected, in which case tan > 52 is preferred. 

In Fig. [4] we show i? yu k versus various sparticle and Higgs masses. In frame a), 
the distribution versus rrih is shown, and we see that low i? yu k prefers the lower range 
of 777/1. This is understandable since low EWFT prefers lower top and bottom squark 
masses, which may not feed a sufficient radiative correction into the computation. 
These low third generation squark mass solutions also tend to give large contributions 
to ^-constraint. If we impose ^-constraints, then rrih can live in the 122 — 128 GeV 
range, at the cost of larger -R yu k- In frame b), we show the distribution versus nig. While 
all solutions — especially low -R yu k ones — favor the heavier range of m g 



2-5 TeV, 

likely beyond LHC reach), the solutions obeying ^-constraints tend to slightly favor lower 
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g, possibly within range of LHC14 searches. In frame c), the distribution versus is 
shown. Here, we see a clear demarcation: 5-constraints favor > 1.5 TeV to suppress 
SUSY loop contributions to BK(B — > X s / -f). This constraint forces the minimum i? yu k 
to move form ~ 1.18 to about 1.3. Additional small flavor- violating contributions to the 
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MSSM Lagrangian could alter the predicted BR(£> — > X s j) and/or BR(£> S — > 
rates and thus allow the lower i? yu k solutions [47]. In frame d), we show the distribution 
versus pseudoscalar. Higgs mass rriA- We see low -R yu k favors the lower range of m^, 
although values up to and beyond 5 TeV are also possible (at the cost of -R yu k ~ 1.4, 
however). Since tan/3 ~ 50, LHC searches for A, H — ^ t + t~ will access a significant 
range of in this case. It is also possible for some range of < 1 TeV for LHC to 
access bA — > bji + n~ production [45] . 

In Fig. [5j we show scatter plots of YUNS points with i? yu k < 1.4 and A < 100 in a) 
vs. mg and b) m~ tx vs. mg space. From frame a), we see as expected that and mg 
are correlated, with some solutions reaching well below m- tx ~ 500 GeV, which should be 
accessible to LHC searches. However, in this case these points all violate 5-constraints, 
so that requiring ^-constraints within measured range requires instead m^, mg > 1.5 
TeV, likely beyond the 14 TeV LHC reach. In frame b), we see also that mg 2 and are 
correlated. This is different from usual NS, where m~ b2 can be far above rni 12 and mg . 
The reason here is that fb is large and so there is also a non-negligible contribution to E u 
from 6 12 . Imposing ^-constraints, we find and mg 2 both > 2 TeV (the former aids 
in lifting into its measured range). 



4 Benchmark points 

In this section we present several YUNS benchmark points, see Table [TJ and compare to 
one YUS benchmark point (HSb, from the "just-so" HS model) from Ref. |43] . 

For HSb, .R y uk ~ 1.02, which is nearly perfect Yukawa coupling unification. The Higgs 
mass rrih — 127.8 GeV is also sufficiently heavy. Unfortunately, the point is now excluded 
by LHC7 searches for multi-jet+i?™ lss plus one 6-tag searches, since rrtg is only 351 GeV. 
We also list in Table [TJ the EWFT measure; for HSb we have A = 2489, indicating 
exceptionally high level of fine-tuning. Much of this comes from the [i parameter which 
turns out to be of order 3 TeV, and thus requires a large value of m 2 Hu at the weak scale 
to cancel against. 

In contrast, point YUNS1 in column 3 has low fine-tuning of A = 39, at the cost of 
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147.4 


112.0 


176.2 


m z, 


49.2 


118.9 


106.2 


170.3 


m h 


127.8 


123.1 


123.8 


123.3 


BR(B X sl ) 


3.1 x 10" 4 


2.7 x 10" 4 


2.8 x 10" 4 


2.4 x 10" 4 


BR(B a -> 


8.1 x 10~ 9 


4.5 x 10~ 9 


4.6 x 10~ 9 


2.4 x 10~ 8 


4613 


0.01 


0.004 


0.007 


c7 SI (Zip) pb 


2.2 x 10~ 13 


4.5 x 10~ 8 


4.7 x 10~ 9 


2.9 x 10~ 9 


a SD (Z lP ) pb 


1.2 x 10~ 9 


7.3 x 10~ 4 


3.6 x 10~ 5 


1.6 x 10~ 5 


{crv)\ v ^o cm 3 /s 


1.9 x 10~ 32 


2.4 x 10~ 25 


3.3 x KT 25 


2.8 x 10~ 25 



Table 1: Parameters and masses in GeV units for HSb [43] and three Yukawa- unified 
natural SUSY (YUNS) benchmark points. We also show 5-decay constraints and dark 
matter relic density and (in) direct detection cross sections. 
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Figure 6: Evolution of Yukawa couplings for benchmark points YUNS1 (left) and YUNS3 
(right) versus renormalization scale Q. 



relaxing -R yu k to 1.37. For YUNS1, rrig ~ 1.3 TeV, with first/secnd generation squarks at 
~ 19 TeV. The point is likely beyond LHC8 reach, but should be accessible to LHC14 
with 10-100 fb -1 . In column 4, we list YUNS2 with -R yu k = 1.3, as low as allowed by 
-B-constraints, but with A ~ 100. This point has m~ g ~ 4.3 TeV and rriq ~ 17.5 TeV, so 
it is likely beyond LHC reach, including a high-luminosity upgrade. While the higgsino- 
like chargino is only 112.5 GeV, it decays via 3-body mode into a higgsino-like Z\ with 
= 106.2 GeV, so that visible decay products are very soft, and likely impossible to 
observe above SM backgrounds at the LHC. Point YUNS3, listed in column 5, features 
R yuk as low as 1.19, with A = 48. This point has BR(B X s j) = 2.4 x 10 -4 , and 
BR(B S — > /i + /U~) = 2.4 x 10~ 8 , so it falls out of the -B-physics allowed range. While gluinos 
and first/second generation squarks are beyond LHC reach, the rather light top and 
bottom squarks may be accessible to LHC searches. All these points have Qzji 2 ^ 0.11, 
leaving room for non-thermal higgsino production and axions. This contrasts point HSb, 
which has a much too thermal abundance and so would need an extremely light axino or 
huge late-time entropy production to tame this over-abundance jUJ. 

In Fig. [6j we show the Yukawa coupling evolution of ft, fb and f T versus renormal- 
ization group scale Q, from m wea k to Mgut for benchmark points YUNS1 and YUNS3. 
These can be compared to similar plots for YUS, as in e.g. Fig. 6 of Ref. [44J. The SUSY 
threshold corrections implemented at the scale Q = m ( - 2 show up as jumps in the 
curves. In the case of YUS models, the mj threshold correction is positive due to a large 
UWj loop contribution which goes like 5b ~ (f^ /327i 2 )([j l A t /m~) tan/3, where both [i and 
A t are extremely large. For YUNS models, with rather low /i, these loops are suppressed 
and injhe case of YUNS1, the gbi loops actually dominate, and are of opposite sign to 
the tiWj loops, leading to the slight downward jump of fb and thus bad Yukawa coupling 
unification. For YUNS3, the UWj loops are larger, and the jump goes upwards, thus 
providing better Yukawa coupling unification. 
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5 Yukawa- unified natural SUSY: LHC, ILC and DM 

searches 



In this Section, we discuss the observable consequences of YUNS for LHC, ILC and direct 
and indirect WIMP and also axion searches. 

5.1 YUNS at LHC 

In natural SUSY models, it is favorable to have multi-TeV first / second generation squarks 
and sleptons because 1. they are safely beyond current LHC searches, 2. they provide at 
least a partial solution to the SUSY flavor and CP problems and 3. they provide additional 
suppression of third generation scalar masses via large 2-loop RGE effects. However, this 
means they are likely beyond any conceivable LHC reach. Third generation squarks may 
be much lighter, and in generic NS models are naively expected to be below the TeV scale 
(but see Ref. [M] where 1-4 TeV third generation squarks work just fine, and lift the value 
of rrih into its measured range). In the case of YUNS, with tan/3 ~ 50, the combined 
light squarks and large tan /3 usually imply violation of ^-constraints, and if these are 
imposed, then top and bottom squarks are beyond 1.5 TeV, and likely inaccessible to LHC 
searches. If additional sources of flavor violation are invoked, then the S-constraints may 
be invalid, and then the solutions with much lower -R yu k ~ 1.2 are accessible, along with 
much lighter top and bottom squarks, potentially accessible to LHC searches. 

For YUNS, the gluino mass may lie anywhere in the 1-5 TeV range. It has been 
estimated in Ref. [45] that LHC 14 with 100 fb _1 should be able to access gluino pair 
production in the case of heavy squarks for rrig up to 1.8 TeV. In the case of YUNS 
with heavier top and bottom squarks, the g is expected to dominantly decay via 3-body 
modes into tiZi and tbWi final states. The gluino pair events will thus contain multi-jets 
plus missing energy plus isolated leptons plus several identifiable fe-jets (46). While the 
higgsino-like chargino and neutralino production cross sections can be large, their decays 
to soft visible particles, arising from the small energy release in their 3-body decays, will 
be difficult to detect at LHC above SM backgrounds [35J. 

In the case where itla < 1 TeV, then it may be possible to detect A, H t + t~, 
especially if these are produced in association with 6-jets, e.g. pp — > bA, bH production. 
It may also be possible to detect bA, bH production with A, H [IB], since the 

production and decay are enhanced at large tan /3. In this case, the A, H mass and width 
may be determined by reconstructing m(/z + yU~). At large tan /3, this width is typically in 
the tends of GeV range and is very sensitive to tan/3. This reaction offers a method to 
distinguish YUNS from NS, in that the former is expected to occur at tan/3 ~ 50. 

5.2 YUNS at ILC 

A linear e + e~ collider operating at ^/s ~ 250 — 500 GeV would in many ways be an 
optimal discovery machine for YUNS. The reason is that by construction \x < 250 GeV, 
so chargino and neutralino pair production should always be available. While the small 
energy release in W\ and Z 2 decay is problematic at LHC, it should be much more easily 
observable in the clean environment of an e + e~ collider. In this sense, an e + e _ collider 
operating at y/s ~ 250 — 500 GeV would be a higgsino in addition to a Higgs factory. It is 
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Figure 7: Thermal neutralino relic density Qh 2 versus Zi mass (left) and versus Z\ higgsino 
fraction (right) for scan points with i? yu k < 1.4 and A < 100. Same color code as in Fig. [TJ 



also possible that some lighter third generation squarks are accessible to ILC with y/s ~ 1 
TeV or CLIC with y/s = 3 TeV, depending if one avoids B-constraints and accepts the 
low mass, low i? yu k solutions. 

5.3 Higgsino-like WIMPs 

A generic prediction of both NS and YUNS models is that the LSP is a higgsino-like 
WIMP with a typical under-abundance of thermally produced (TP) neutralinos f2~ p /i 2 ~ 
0.002 — 0.01. However, in cases where m\/2 is as low as ~ 300 GeV and [i is as large as 
200 — 250 GeV, then there can be substantial bino-higgsino mixing, boosting Q]~ p h 2 up 

where we plot rf~ p h 2 versus 

nig and versus the Z\ higgsino fraction from YUNS models satisfying i? yu k < 1.4 and 
A < 100. 

The typical under-abundance is an appealing feature if one invokes the Peccei-Quinn 
solution to the strong CP problem, in which case one must introduce an axion superfield 
a which contains a pseudoscalar axion a as well as a spin- 1/2 axino a and a spin-0 
saxion s. In such models, one expects both saxion and axino masses at or around the 
SUSY breaking scale, so that dark matter is comprised of an axion- WIMP admixture. 
In this case, thermal production of axinos and subsequent decay to states such a gg 
bolster the WIMP abundance beyond its TP-value. In addition, axions are produced as 
usual via coherent oscillations. It is also possible to suppress the WIMP abundance in 
the mixed aZ\ cosmology via late-time entropy production from saxion production and 
decay, although this case tends to be highly constrained by maintaining successful Big 
Bang nucleosynthesis (BBN). The upshot is that in the mixed aZ% dark matter scenario, 
it may be possible to detect both a WIMP and an axion. 

In the case of the YUNS model, the higgsino-like neutralinos have a substantial spin- 
independent (SI) direct detection cross section, as illustrated in Table [l where a sl (Zip) ~ 
10 -8 pb. While this level of direct detection cross section is now highly constrained by 
recent XENON100 [49J results, one must bear in mind that the higgsino-like WIMPs 
would constitute only a portion of the dark matter, so their local abundance might be 
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up to a factor of 30 lower than is commonly assumed. In Table [lj we also list <r SD (Zip) 
(relevant for WIMP detection at IceCube) and (av)\ v ^o, relevant for detection of dark 
matter annihilation into gamma rays or anti-matter throughout the cosmos. While these 
cross sections are also at potentially observable levels, again one must take into account 
that the overall WIMP abundance may be up to a factor of about 30 below what is 
commonly assumed. Plots of a sl (Zip), a SD (Zip) and (crv)\ v ^ have been presented in 
the case of higgsino-like WIMPs in Ref 's [35j and (37] and so similar plots will not be 
reproduced here. 

6 Summary and conclusions 

Previous analyses of t — b — t Yukawa-unified models suffer from two problems: 1. they 
tend to predict a light gluino rrig < 500 GeV (alhough solutions are possible for much 
heavier gluinos) which is now excluded by LHC searches, and 2. they suffer from extreme 
fine-tuning in the electroweak sector. In this paper, we examined how well the Yukawa 
couplings could unify in the natural SUSY context, where /i ~ 100 — 250 GeV, while at 
the same time requiring the light Higgs mass trih ~ 122 — 128 GeV. The small value of // 
suppresses the large UWj loop contributions to the 6-quark Yukawa coupling which seem 
to be needed for precision Yukawa coupling unification. Nonetheless, by scanning over 
NUHM2 parameters with split third generation, we are able to find solutions with -R yu k as 
low as 1.18. These solutions, with very light third generations squarks and tan /3 ~ 50 tend 
to violate I?-physics constraints. If 5-physics constraints are imposed, then only R ~ 1.3 
can be achieved. However, the I?-physics calculations can be modified if additional small 
flavor- violating terms are allowed in the MSSM Lagrangian, so it is not clear how seriously 
R > 1.3 should be taken. 

The Yukawa-unified natural SUSY spectra have important differences from previous 
YU spectra. The gluino mass can easily be in the 1-4 TeV range, thus avoiding LHC 
constraints from SUSY searches. The light higgsino-like charginos and neutralinos decay 
to soft particles, also avoiding LHC searches. However, light higgsinos should be easily 
accessible to an ILC with yfs ~ 0.25 — 05 TeV, as is typical of all NS models. In 
addition, as in all NS models, the lightest neutralino is higgsino-like with a typical thermal 
underabundance of WIMP dark matter. We regard this as a positive feature in that the 
WIMP abundance is typically increased in non-standard (but more attractive) cosmologies 
such as those conatining mixed axion-neutralino cold dark matter. 

The question arises as to how to distinguish YUNS from ordinary NS. The YUNS 
model requires tan ~ 50, which leads to large production cross sections for heavy 
Higgs bosons A and H at LHC, and large widths for these particles. If the rare decays 
A, H fi + fi~ can be identified with suficiently high statistics (perhaps at a luminosity 
upgraded LHC), then the widths may be measured with precision, allowing one to highly 
constrain tan (3, and perhaps verify that it is consistent with YUNS models. 
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